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LAMINAR SLIP PLOW IN A FLAT DUCT OR A ROUND 


TUBE WITH UNIFORM WALL HEAT TRANSFER 
"by Robert M. liman 
Lewis Research Center 

National Aeronautics and Space Administration 
Cleveland, Ohio 

ABSTRACT 

An analysis was made to determine the effects of low-density phe- 
nomena on the heat-transfer characteristics for laminar flow in a 
parallel-plate channel (flat duct) or in a circular tube with uniform 
wall heat flux. Consideration was given to the slip-flow regime, 
wherein the major rarefaction effects are displayed as velocity and 
temperature jumps at the conduit walls. The results obtained apply 
along the entire length of the conduit, that is, in the thermal entrance 
region as well as far downstream. The solutions contain a series ex- 
pansion, and analytical expressions for the complete set of eigenvalues 
and eigenconstants for this problem are presented. The results give 
the wall temperatures, Nusselt numbers, and thermal entrance lengths for 
the conduits for various values of the rarefaction parameters. 
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NOMENCLATURE 

a accommodation coefficient 

coefficient in series for temperature distribution in parallel- 
plate channel 

coefficient in series for temperature distribution in circular 
tube 

specific heat at constant pressure 


m 


2 


D n 

Dj 

d 

f(n) 

f(®) 

Q(ti) 

g 

H(a>) 


II, 1 ! 


L 

2 

M 

H 

Hu 

Pr 

P 

q. 

R 


coefficient defined lay equation (32) 
coefficient defined "by equation (57) 
thermal diameter, 8L/c 
tube diameter, 2rQ 
constant defined in equation (31) 

dimensionless velocity for parallel-plate channel, u(tj)/u 
d ime nsionless velocity for circular tube, u(co)/u 
transverse temperature distribution in fully developed region 
for parallel -plate channel 
specular reflection coefficient 

transverse temperature distribution in fully developed region 
for circular tube 

heat-transfer coefficient, qj (ty- - tfc) 

value of definite integral, equation (34) for par all el-plate 
channel, equation (59) for round tube 
Bessel function of first kind and first order 
half distance between plates 
mean free path 

constant defined by equation (60) 
constant defined by equation (61) 

Husselt number, hDrp/je or hd / K 
Prandtl number, \iC^/k 
static pressure 

rate of heat flux per unit area from wall to fluid 
transverse or radial distribution function 
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Re Reynolds number, 2p?CL/u for parallel-plate charnel, pTO/p 

for circular tribe 
Rg gas constant 

eigenfunctions of equation (14) for parallel-plate channel 
% eigenfunctions of equation (43) for circular tribe 

r radial coordinate 

Tq tribe radius 

t temperature 

tg gas temperature adjacent to wall 

u reloelty 

x axial coordinate 

y transrerse coordinate 

Greek symbols: 

a dimensionless Telocity slip coefficient, i u /2L or | u /d 

ft- ^ /2l i 


m 

% 




T ratio of specific heats 

£ dimensionless coordinate, x/2L or x/r^ 

t] dimensionless coordinate, oy/2L 

k gas thermal conductlrtty 

\ separation constant 

eigenvalues of equation (14) for parallel-plate channel 
7 % eigenvalues of equation (43) for round tribe 

p absolute Tlseoslty 

^ temperature-jump coefficient 

yeloclty-sllp coefficient 
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p gas density 

er symmetry ntnriber 

<p rarefaction parameter, p or p ^/Rgt/pd 

i|f dimensionless quantity, RePr/ 0 2 for parallel-plate channel, 

ReRr/4 for circular tube 

03 dlmensionlesB coordinate, r/r^ 

Subscripts] 

b bull condition of gas 

d fully developed region 

d,c fully developed region for eoutinuum flow 

e entrance region 

1 gas entering channel, x «= 0 

s slip condition at vail 

v vail 

0 heated section entrance, x = 0 

Superscript* 

*™ average value 

INTROEUCTION 

In recent years, considerable interest has developed in the study 
of the fluid-flow and heat-transfer characteristics of rarefied gases* 

This interest has been stimulated by the Increasing frequency of low-dens lty~ 
environment applications and the advent of space flight* Only very re- 
cently has attention been directed to the problem of heat transfer to 
rarefied gas flow in conduits* 

Of particular Interest in internal, rarefied gas-flow studies has been 
the problem of lanrtun-r forced-cotrrection heat transfer in conduits under 
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slip-flow conditions [l and 2 ]* The essent i al simplifications introduced 
In these investigations to obtain analytical eolations are fully estab- 
lished temperature profiles ana fully developed velocity distributions* 

The present investigation is concerned -with, the more general problem 
of determining the heat-transfer characteristics along the entire length 
of 'tiie conduit, that is, in the thermal entrance region as veil as far 
downstream, for laminar slip flow in a parallel-plate channel (often re- 
ferred to as a "flat dnet") or In a circular take with uniform wall heat 
flux* 

In the section FLOW IN EAEAIIiEL^HATE CHANNEL is considered the prob- 
lem of slip flow of a rarefied gas In a parallel-plate channel with uni- 
form -wall heat flux at one or at both walls* Both heating arrangements 
are frequently encountered in practical applications. The problem of slip 
flow in a circular tube with uniform wall heat flux is discussed in the 
section FLOW IN CIRCUIAR TUBE# 

FLOW IN PARATJiFL-ELATE CHANNEL 

The coordinate systems for the problems under study are shown in 
Fig, 1# A slightly rarefied gas flows in the positive x-dlrectlon with 
a fully established velocity profile# Up to a point (x * 0) the channel 
walls and gas are isothermal at temperature tj_* After this point a uni- 
form wall heat flux is applied* It is desired to determine the tempera- 
ture distribution and the variation in the heat-transfer coefficient along 
the entire length of the channel* 

It is convenient to place the plane y = 0 at the plane of symmetry, 
that is, at the middle of the duct in the case of heating of both walls at 


y = ±L (Fig* 1(a)) and at the insulated wall in the ease of heating from 
cote side at y = +2L (Fig* 1(h))* Beth eases are included in the follow- 
ing development by defining a symmetry number o', which is also the 
nuniber of heating surfaces [3]* 

Before the energy equation can be solved, the gas Telocity distri- 
bution mast be known* This distribution -was investigated in [4]* The 
use of the results leads to the dimensionless velocity profiles u( t]) fit 
as 

f ( T)) = ( 3 / 2 ) (1 - T)^ + 4a) /(l + 9fc) 0 m 2 (l) 

f ( T|) = 6(r) - II 2 + »)/(! + 6tt) 0 = 1 ( 2 ) 

where a a f^/21% The slip coefficient 1 b given by the expression 
[5] 


&tt “ t(2 - g)/g]Z (3) 

where Z is the mean free path and g is the specular reflection co- 
efficient* The relation between the average velocity and the slip ve- 
locity is easily obtained as 

u^yfa ** 6a/(l + 6a) o = 1,2 (4) 

The starting point of the heat-transfer analysis is the differential 
equation for convective heat transfer in the parallel-plate channel flow 
with fully established velocity profile# With the gas properties 
assumed constant, the heat conduction in the flow direction compared 
with that in the transverse y-dlrection assumed negligible, and the 
viscous dissipation assumed negligible, the equation can be written in 
the form 


pC u(dt/cfcc) = ^(cJ^t/cy 2 ) 


(5) 
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Eq^atloai (5)* written In terms of dimensionless -variables, ■becomes 

*♦ (^/bii 2 ) (6) 

The boundary conditions are as follows* 

Uniform wall heat flux* 

St/Si] *= (2L/s)(q/fc) at ii ■» 1# t, J> 0 

Symmetry* 

dt/^rj » 0 at t] » 0 
Specified entrance tempera-fore* 

t at £ = 0 

Whan the wall heat flux 1 b uniform* It la known that for Tory large 
-mines of x there 1 b a fully developed thermal situation characterized 
hy a linear rise In the temperature at all point* In the cross section 
along the channel^ that is* 

(dt^/chc) = (oq/RePr/e) a * 2*1 (7) 

Equation (7) can he written alternatively as 

(t a *- t 1 )/(2I«/<r)(q//c) = O 2 ( x/ 2ti) /ReEr + G( iq) (8) 

The function G( t]) for each value of C Is 

G(ti) - [(3/4) ri 2 - (1/8) ti 4 - (39/280)1 + [-(l/4) T] 2 + (l/8)ri 4 

- (13/280) Kttg/t) + [2/1051(^^0 2 ff - 2 (9a) 

G(ri) « [q 3 - (1/2) T] 4 - (9/70)1 + [(X/2) iq 4 - r\ Z + (1/2)t) 2 

- (s/rontug/u) + [1/210I (ug/u ) 2 <r« i ( 9 b) 

The details of the analysis for G(t]) are omitted In the present study 
hut may he found, in [4]* The first quantity In brackets on the right side 
of (9a) and (9b) represents -foe usual transverse temperature distribution 
for continuum flow conditions* while the second and third quantities in 
brackets are connected with one effect of gas rarefaction, namely* that 
of velocity jump* 
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Determining the soLLutlcm in. the entrance region 1* cenrenient If 
a difference temperature t e Is defined as 

y 5*t|) - t(CfTi) (10) 

23isn t e meat satisfy the relation 

W(ri)(te 9 M) - (31) 

•with boundary conditions 

btg/^T) *» 0 at t) = 0 and tj = 1 (12) 

The solution of (11) that -will satisfy (12) can "be found hy using 
a joroduct solution that leads to a separation of TBrlablee* This solu- 
tion will hare the fora 

t e /(2L/o)(q Jk) » Y\ C m R m ( T]) Q3Cp[*-o 2 \ 1 (x/2L)/RePr] (13) 


where 7^ and are, respectlYely, the elgemralues and eigen- 
function* of the Sturm-Liorccvllle problem j 

(d 2 R/dri 2 ) + 7f(ti)R « 0 
dR/di) b 0 at t) ■ 0 and r\ » 1 
The coeffielents C m In (13) are determined so as te satisfy the 
"boundary condition at the entrance to the heated section £ » Ot 


a * 2,1 


(14) 


t(0,T])/(2l/ff)(qyK) - -G(ti) 




(15) 


This result together with the orthogonality property of the eigenfunctions 
leads to 




^ G(rj)f( T]) R m ( ti) dii f(ti)'lg(Ti)dTi 


J' 1 G(ti)f(r 1 )R ffi (Ti)dr 1 / [RU)d 2 R/^ * * 2 A 

0 ( 16 ) 
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The Integral la (16) can he eraluated hy substituting g(t]) from (9a) on 
(Sb) and f(T)) from (1) or ( 2 )* Integrating "by parte,’ and utilizing (14) * 
The final result Is Identical for hath cases j 

Cm “ 1]=!^)^=)^ 0 “ 2 A (17) 

• The. complete seintion for the tenrgera t are that applies in hath the 

entrance and the fully denreleped regions is found by adding the solutions 
for t fl and t 0 * 

A result Of practical interest is the longitudinal Tarlatlon of wall 
temperature t^ corresponding to a uniform wall heat flux# Before the 
wall temperature warlatlon can he determined, howerer, it is necessary 
to consider another effect of gas rarefaction that enters through the 
thermal boundary condition at the wall; which permits a Jump between the 
surface temperature t^ and the adjacent gas te m p erat ure t g [5] j 

"M^ArW/® (18) 

where represents a temper a t ur e-jump coefficient related to other 
properties of the system hy 

it « 1(2 - a)/a][2r/(r + 1) K *A*) (is) 

For uniform wail heat flux, 

&M^ZL/o “ E(^t a /^y) + (&t*Ar) 1^21/0 - (flA) 

so that the tem pera t u re Jump at the wall ran he written 

t g - t v * -(2 L/o) (q/ff) [s( |^./2L) ] (20) 

Ccmhlnl ng (8) and (13) in accordance with (10), setting tj = 1, and 
then combining the result with (20) yield 
(\ - t 1 )/(2L/cr)(q//c) - C 2 (x/2L)/RePr + G(l) 

00 

+ 0(^/2!) + 23 C m R m ( 1) exp[-c2y x/2L) /RePr ] (21) 
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For a uniform wall heat flux, the hulk t em p er ature t^x) 1 b given "by 
% - t± + ff 2 [ ( 2L/ff) ( oJk) (x/2L) /RePr] (22) 

Combining (21) and (22) yields 
(\ - t b )/(2L/o)(g/K) = G(l) + «r(^/2L) 

21=1 

Then, for the folly developed situation (x -* <»), 

(\ ** t b ) a/ /(2L/ff)(^) » G(l) + 0(^/2!) (24) 

Dividing (23) by (24) yields the Important ratio 

(t * ‘ V/ft » ' - {° (1) + 0(H/21) 

+ £ C ffi I^(l)e X p[- C r 2 ^ n (x/2L)/RePr]l/[G(l) + *0^21)] (25) 

m=l J 

Equation (25) out be evaluated once numerleal values of 7^, 1^(1), and 
Cpj bare been obtained for given values of a or Ug/Eu 

The Nusselt number may be determined from the definition of a thermal 
d i am e ter Drp, which depends on the area of the heating surface [3]* Far 
the present analysis. Dtp a 8L/c* Then, when this definition Is applied, 

3m = iuji^/k = [q/(ty ** tj) ) ] ( Sl/o/C) 

Wlaen (23) Is used. 


(23) 


Hu = 4 /|g(1) + 0(^21) + c m \ ( 1) X / 2L ) /RePr 1 j> (26) 


It Is of interest to examine the belaavior of the Nusselt number 


at the entrance of the heated section (x = 0) and also in the fully 
developed ration (x ~). At the entrance. 
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G(l) + o(l t /2L) + YZ 0^(1) 
m«l 


From (15), however, setting n = 1 yields 



n «0 ■ 4/cr(I t /2L) cr = 2,1 (27) 

In tile absence of a temperature- Jump effect, the local Nusselt number 
starts with NUq -*• °°» With a temperature jump, however, the local Husselt 
number commences with a finite value given by (27)* The effect of the 


number of heating surfaces enter b through the symmetry number a* 

When (26) is used, the fully developed Nusselt number is deter- 


mined as 


- 4/[G(l) + 0(^/2 L)] a « 2,1 (28) 

The fully developed Nasselt number becomes for each case 
No* - (140/17 )/[l - (6/17) (ttg/u) 

+ (2/51) (Ug/u) 1 2 + (70/17) (lt/2L)] a = 2 (29a) 

WOd * (140/13)/[1 - (3/26) (u g /u) 

+ (1/78) (Ug/u ) 2 + (35/13) (I t /2L) a = 1 (29b) 

In the absence of rarefaction effects, the fully developed Husselt number 
has the value 140/17 .« 8,23 (a = 2) or 140/13 = 10,77 (a « l). 
From (29a) or (29b), it is clear that the effects of velocity jump (Ug / 0) 
would tend to increase the Nusselt number, while the temperature jump 
would act to decrease the Husselt number. S imilar results have been ob- 


served in the case of circular-tube slip flow [1,2], 

numerical values of the entrance Nusselt number (27) have been 
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evaluated at functions of the parameter uV^gV^P 1, and. related to the 
mean free path l by the relation [5l.pyR g t/2pL « *^/2/x (z/2L), The 
raises are plotted in Fig* 2* The raises for r and Prandtl number are 
representative of air and moat diatomie gaaea* Fully dorr-eloped Nusselt 
mm* era aa glren by (29) hare been erraluated aa a function of the param- 
eter ji-y'Rgt/ 2pL and are plotted in Fig* 3 for X ** 1 * 4 , Pr * 0 , 73 , 
g * If and a * 1*0 and 0,4 in the form of the ratio Nu,j/Ku^ c 
■where nu^ 0 ia the appropriate fully developed continuum raise, The 
effect of gaa rarefaction ia almya to decrease the ralue of the Nusselt 
number below its continuum ralue, the result being more pronounced for 
two-aided heating than for one-sided heating* 

The d imen sionless wall-to-bulk temperature difference ’d 


(ty - " tfc) , (25), and the local Nusselt number No, (26), can 

he erraluated along the entire duct length aa aeon aa the eigenvalues 
Aju, eigenfunctions R m ('n), and series coefficients ^ hare been deter- 
mined, The function R^ti) is the solution of (14) t 

(d 2 R/dn 2 ) + Xf(ii)R - 0 

( 

dR/dii “0 at •»! « 0 and i) « 1 
The normalization conrention 


2,1 


(30a) 


R(0) - 1 


(30b) 


is also used* 

Asymptotic expressions for the eigenvalues 7^ and constants 
and I^(l) for both the case o « 2 and the case a ** 1 are derired in 
[4] and the results are presented aa follows! 

tan » | IfZ + (1 + 4a)sin“ 1 Cl/-yi 7 + 4a)] /2a (4a) 3 / 2 = E a (31) 
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\ 3 <VW(1) - *8«/[e a + 1 + Oj/E 0 )] (32) 

vhsrc 

P m S *V^ Ii (33) 

Vf(n) ^ - V3/8 [i/C + (1 + 4a)ain“' 1 (l/Vl + 4*)] jw*- + 6a 

(34) 

and E 0 , fear a given, value of % la a constant* Equation (31) indicates 

that E^ * 2Eg* Tabulated values of the first fire roots of (31) for a 

number of values of E 0 are given In [6]* The values of 1^ for any 

given slip velocity ratio Ug/u are shown in Fig* 4* 

The first four eigenvalues and coefficients are shown In table l(a) 

for symmetrical two-sided heating (a = 2) or In table 1(b) for unsym- 

metrical one-sided heating (o ■ l). The results for continuum flow were 

obtained from expressions presented In [3], while for slug flow) the 

eigenvalues were obtained as the positive roots of sin 7^/^ = 0 or 
l/2 

» ms. The coefficients were obtained from the equally simple 

result % » -2/}^, 

The numerical value of l/3 for Ug/u corresponds to ^/2L * 0.0833, 
while a value of 0/5 for Ug/u corresponds to l^/EL = 0.25. The results 
for Ug/u = 1 (slug flow) are outside the slip regime but have been in- 
cluded as limiting values and for comparison. 

The level of accuracy of the foregoing results was checked [4] by 
computing the eigenvalues and eigenfunctions of (30a), as well as the 
coefficients given by (16), on an electronic (IBM 7094) computer, 

by the Runge-Kutta method, for Ug/u = l/3 and 3/5. The eigenvalues and 
coefficients so obtained are listed in table I. The relevant quantities 
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as computed from the previously presented analytical expressions are in 
remarkably close agreement with the electronically computed -values, es- 
pecially for values of m ^ 2. It is concluded that the asymptotic 
formulas are suitable for m ^ 2, 

The variation of the dimensionless wall-to-bulk temperature difference 
along the duets can he evaluated with the numerical information given in 
table I, Before proceeding with the evaluation, however, it is illuminat- 
ing to examine the wall-to-bulk temperature difference at the entrance of 
the heated section. With x = 0, (25) becomes 

- *b> a - [°(i> + 

00 

m**l 

When H m 1* hfirwenrer* (15) becomes 

G(l) « - ^ N / ^1^(1) 
m*l 

so that 

(t v - t- b ) Q j(t v - t-b) d - a(l t /2L)^|G(l) + c0t/2L)] (36) 

Tn the absence of a teo^e:rature*jump effect* 'bite mll*to*bulk tempera* 
tare difference is zero at the entrance* With, a temperature jump* how* 
ever* the entrance temperature difference has a nonzero value* Equation 
(36) is plotted in Big. 5 as a function of the two parameters u^/u and 
tfc/2L for a = 2 and 1» For either wall heating situation the entrance 
temperature difference increases with an increasing value of &^/2L, 

The magnitude of th© slip velocity has only a small influence on the 
quantity (t w - - tb) d for 0 - 1/ while for a - 2 the influence 

of slip velocity is more pronounced. 


[g(1) + o(!t/2L)J (35) 
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The Tar 1st Ion of the wall-to-bulk temperature difference along the 
duct length ms evaluated from (25) for several values of the rarefac- 
tion parameters %/u and 2L» Plots are given in Figs. 6 and 7 

for a » 2 and a » 1* respectively. 

Inspection of Figs. 6 and 7 reveals several Interesting trends. 

First of all, for a fixed value of i^f 2L, the mil temperature varia- 
tion Is more sensitive to the slip velocity over most of the duct length 
for the unsymmetrically heated channel (a = l) than for the symmetrically 
heated channel (a = 2). Near the entrance, however, the reverse effect 
is obtained. For both mil heating situations, the slip velocity has the 
effect of retarding t w - t^ in its approach to the fully developed 
value, while the temperature jump has the opposite effect. Finally, the 
abscissa scale for a » 1 Is twice that for a = 2, The length required 
for the wall-to-bulk temperature difference to approach fully developed 
conditions is thus greater for the unsymmetrically heated c hann el than 
for the channel heated uniformly from both mils. 

It is the practice to define a thermal entrance length as the heated 
length required for t w - t-^ to approach within 5 percent of the full y 
developed value. A horizontal dashed line corresponding to an ordinate 
of 0.95 is shown in Figs. 6 and 7. 

It is perhaps somewhat more i lluminat ing to present the variation 
of the wall-to-bulk temperature difference in terms of the rarefaction 
parameter |i-/R g t/2pL as in Fig. 8. The effect of increased gas rare- 
faction is to shorten the thermal entrance length. The accommodation 
coefficient also has an important effect on the thermal entrance length, 
and this effect is associated with the increase in temperature jump with 
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decreased accommodation coefficient. 

Tile longitudinal variation of the Nusselt number along the duct 
with uniform heat flux at one or both mils (26) ms evaluated and is 
plotted in Fig. 9. The velocity and temperature jumps give rise to 
opposite changes in their effect on the WusBelt number variation* the 
velocity Jump tends to increase the Nusselt number at a given axial 
position, "while the temperature jump tends to decrease the Nusselt 
number. Numerical evaluations of the Nusselt number dependence on the 
parameters [i-v'Rgt/2pL a ad a are plotted in Fig* 10, Clearly the 
overall effect of the gas rarefaction is almys to decrease the Nusselt 
number below its continuum value at every position along the heated 
length, 

FLOW IN CIRCULAR TUBE 

Attention is now turned to the case of axially symmetric slip flow- 
in a circular tube. The coordinate system for the present problem is 
Shown in Fig. 11, 

It is again assumed that the velocity profile is fully developed 
and is unchanging along the tube length. The velocity distribution and 
the fully developed heat -transfer characteristics have already been in- 
vestigated for the circular-tube case [l], and many results obtained 
are Immediately applicable* The development of the round-tube system 
is similar to that of the parallel-plate channel. 

The differential equation for convective heat transfer is now 

p CpU (dt/dx ) «■ ( fc/r ) (d /dr ) ( r St /dr ) ( 37 ) 

The assumptions and restrictions of this equation are the same as those 
previously explained. When written in terms of the dimensionless vari- 
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ables, (37) "becomes 

2tf (os) (dt/c)£) * (l/oo) (c)/cko) (cd dt/ck») (38) 

The Telocity distribution u is given in [l] and from the restats the 
dimensionless Telocity profile f(co) and the slip Telocity ratio Ug/u 
are easily obtained: 

f(m) = 2(1 - cn 2 + 4a)/ (l + 8a) a a l^/d (39) 

Ug/u = f(l) = 8a/ (1 + 8a) (40 ) 

For large x, a fully developed temperature profile t d exists in the 
form 

(t d - t i )/(q.r 0 /K) = (4/RePr) (x/'r 0 ) + H(cd) (40) 

■where the radial ftmcfcion H (cd) Is given In [ll as 

H(a>) - co 2 - (l/4)co 4 - (7/24) - £i/2)cd 2 - (l/4)co 4 J(u B /u) + (l/24) (ug/u) 2 

(41) 

The solution for the thermal entrance region can be shown to haTe 
the form 


t e / (qr 0 /K) = C n R n (©)exr[-4?s n (x/r 0 )/RePr] (42) 

where Tv^ and Rg^ are, respectively, the eigenvalues and eigenfunctions 
of the Sfcurm-Liouville problem: 

(d/cko) [co(dR/dcn)] + 27tfof(co)R = 0 
cLR/cLd = 0 at co = 0,1 

The coefficients C n in (42) are obtained from the result 


C n = - 


'i , r 1 o 

2coH(co)f (oo)R n (co)dn) / / 2oof (cD)R n (a>)d£D 

J 0 I JO 


2coH(oo)f(co)R ] 


: n (a>)dco/| 


[R(<d) (d Rfia &X)] 


oo=l 


(43) 


(44a) 


n 
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or 


c n ■ 1 ^[>v(S 2 r/^d 

In order to obtain (44a), the result 


(44b) 


H(m) 


00 

- - 52 C n R n(“) 

n«=l 


•was used. 

The complete solution for the temperature that applies along the 
entire tube length is obtained by summing (40 ) and (42) to obtain 
(t - t 1 )/(qr 0 //c) = 4 (x/ r 0 ) / RePr + H(co) 

+ ^2 c n R n (co) exp [ -4?^ (x/r 0 ) / RePr ] (46) 

The temperature-jump effect at the tube wall is given by 

^g - t v = -2(q.r 0 //c)(| t /d) (47) 

Hence, the wall, temperature along the length of the tube is obtained as 
(t w - ti)/(q.ro//c) = 4(x/r 0 )/(RePr) + H(l) + 2(| t /d) 

00 

+ V. C n R n ( 1 )exp[-4X n (x/ro)/ (RePr)l (48) 

n«l 

This equation can be rephrased in terms of the bulk temperature tj,(x) 
with the result 

r 

(tv * tb)/(t v - tb) d - E(l) + 2(| t /d) 

+ ^ C n R n (l)esp[-4X n (x/r 0 )/RePr]\ /[h( 1) + 2 (fit/d)] (49) 


n=J 


where 


and 


(t v - t^) - (qr 0 //f)[H(l) + 2(l t /d)] 

a 


(50) 


t|j(x) = t ± + (qr 0 / k) (4x/r 0 )/RePr 


(51) 
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The' Nusselt number may be determined from the definition 

Nu = h(2r 0 )/< = [ g/(tw - t b )](2r 0 /k) 
When (49) and (50) are used, the result obtained is 


Nu = 2 /|h( 1) + 2(| t /d) + g C n R n (l) eX p [4A n ( X /r 0 )/RePr] 


( 52 ) 


(53) 


The Nusselt numbers at the entrance of the heated section Nuq and in the 
fully developed region Nu d are readily obtained from (53) as 

Nu 0 = l/(g t /d) (54) 


Nu d = (48/ll)/[l - (6/ll) (Ug/ti) + (l/ll)(u s /Tl) 2 + (48/ll)(| t /d)] (55) 

The effects of gas rarefaction on the fully developed Nusselt n umb er (55) 
have been considered in [1]. 


The asymptotic behavior of (43) at large values of ~K is examined 
in [4,7] and it is shown that the asymptotic expressions for the eigen- 
values A n and constants D n = C n R n (l) are given by 


+ N 


wnere 


taJa Mp n - N 
D n = -16a/[N + (N 2 /M) + Mpf] 

^n = ^l 


(56) 

(57) 

(58) 


-1 


= V2fIcuT (to = + (l + 4a)sin _1 (l/- v /l + 4a)]/vi + 8a 


M = 4 



(59) 

-^ia + (1 + 4a)sin" 1 (l/-,yi + 4a )J (60) 

N = 1 - 4a (61) 


The values of I b for any given slip velocity ratio u s /n are shown in 
Fig. 12. 
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The first four eigenvalues and coefficients for the case of flow in 
a round tube are shown in table II. The results for continuum flow 
(u s /n = 0) were obtained from expressions given by Dzung [8], while for 
slug flow (u s /u *♦ l) the eigenvalues were obtained as the roots of 
Jl(2A n ) 1 / 2 = 0. The coefficients D n are then obtained from the simple 
result D n = -l/A n . Also shown in table II are the data obtained through 
the use of an IBM 7094 computer by the Runge-Kutta method [4], It is 
apparent that the asymptotic formulas yield values of sufficient accuracy 
for n > 2. 

Numerical values of the Nusselt number variation along the tube 
length (53) were evaluated as functions of the two parameters u /n and 
l^/d and are plotted in Fig. 13. The trends are similar to those ob- 
served in the parallel-plate channel system. The Nusselt number varia- 
tion 1 can be calculated as a function of the parameter p-/R t/pd; the 
results of such a calculation are plotted in Fig. 14. Increased gas rare- 
faction and/or decreased accommodation coefficient reduces the Nusselt 
number below its continuum value, and, in addition, shortens the thermal 
entrance length, which has been defined alternatively as the heated 
length required for the Nusselt number to approach within 5 percent the 
fully developed value as given by (55). 

OTHER RAREFACTION EFFECTS 

In [4], modification of the present heat-transfer results for laminar 
channel or tube slip flow is made, or discussed, to account for additional 
slip effects, such as wall shear work, modified temperature jump, and 
thermal creep velocity. 
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CONCLUSIONS 

Solutions were obtained for laminar, forced- convection heat transfer 
to a slightly rarefied gas flowing between parallel plates or in a circu- 
lar tube with uniform wall heat flux. The wall temperatures and Nusselt 
numbers in the entrance and fully developed regions can be obtained as 
functions of the velocity and temperature jumps at the wall, or as func- 
tions of the mean free path. 

The results indicate that the slip-flow Nusselt numbers are lower 
than those for continuum flow at all, axial locations along the conduits 
and also that the thermal entrance length is decreased with increased 
gas rarefaction for either the parallel-plate channel or the circular tube. 
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TABLE I. - EIGENVALUES AND COEFFICIENTS FOR SLIP FLOW IN A 
PARALLEL-PLATE CHANNEL WITH UNIFORM HEAT FLUX 
(a) Symmetry number, 2. 




Ratio of slip to average velocity, u s /u 

0 

1/3 

3/5 

1 

Analytical 

solution 

Numerical 

solution, 

Runge-Kutta 

method 

Analytical 

solution 

Numerical 

solution, 

Runge-Kutta 

method 

Eigenvalue 








3.540 

3.78 

3.33 

3.35 

3.23 

3.14 

W 2 

6.800 

6.72 

6.49 

6.41 

6.36 

6.28 


10.05 

9.78 

9.65 

9.54 

9.50 

9.42 

-k\/ 2 

13.30 

12.90 

12.82 

12.69 

12.65 

12.56 

Coefficient 







E 1 

-0.2090 

-0.1479 

-0.2331 

-0.2110 

-0.2264 

-0.2030 

D 2 

- .0703 

- .0642 

- .0701 

- .0613 

- .0618 

- .0508 

d 3 

- .0367 

- .0332 

- .0336 

- .0282 

- .0280 

- .0226 

D 4 

- .0230 

- .0198 

- .0197 

- .0165 

- .0161 

- .0127 


(b) Symmetry number, 1. 




Ratio of slip to average velocity, u s /u 

0 

1/3 

3/5 

1 

Analytical 

solution 

Numerical 

solution, 

Runge-Kutta 

method 

Analytical 

solution 

Numerical 

solution, 

Runge-Kutta 

method 

Eigenvalue 








3.800 

4.09 

3.50 

3.51 

3.35 

3.14 


7.071 

6.99 

6.66 

6.51 

6.46 

6.28 

^ 2 

10.33 

10.01 

9.82 

9.61 

9.58 

9.42 


13.60 

13.09 

12.98 

12.72 

12.71 

12.56 

Coefficient 







D 1 

-0.1470 

-0.0711 

-0.1821 

-0.1685 

-0.1920 

-0.2030 

»2 

- .0525 

- .0425 

- .0583 

- .0567 

- .0566 

- .0508 

D 3 

- .0278 

- .0259 

- .0291 

- .0271 

- .0267 

- .0226 

D 4 

- .0176 

- .0169 

- .0175 

- .0156 

- .0154 

- .0127 















TABLE II. - EIGENVALUES AND COEFFICIENTS FOR SLIP FLOW IN A 


CIRCULAR TUBE WITH UNIFORM WALL HEAT FLUX 




Ratio of slip to average velocity, u s /u 



0 

2/5 

2/3 

1 



Analytical 

solution 

Numerical 

solution, 

Runge-Kutta 

method 

Analytical 

solution 

Numerical 

solution, 

Runge-Kutta 


Eigenvalue 







2.531 


2.55 

2.64 

2.60 

2.710 


4.578 

4.71 

4.63 

4.75 

4.74 

4.955 

* \ ,z 

6.599 

6.76 

6.69 

/. 6.88 

6.86 

7.195 


8.610 

8.81 

8.75 

9.00 

8.98 

9.425 

Coefficient 







D 1 

-0.1985 


-0.1855 

-0.1670 

-0.1658 

-0.1360 

D 2 

- .0693 

-0.0594 

- .0605 

- .0515 

- .0510 

- .0406 

D 3 

- .0365 

- .0306 

- .0301 

- .0247 

- .0245 

- .0194 

D 4 

- .0230; 

. - .0217 

- .0185 

- .0145 

- .0144 

- .0113 
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(a) Symmetry number, 2. (b) Symmetry number, 1. 

Figure 1. - Physical model and coordinate system for parallel-plate channel. 


Accommodation 

coefficient, 

a 



Figure 2. - Effect of gas rarefaction on entrance Nusselt 
number for parallel -plate channel. 
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Accommodation 

coefficient, 

a 


0.4 Symmetry 

1.0 number, 



Rarefaction parameter, /x^Rgt/2pL 

Figure 3. - Fully developed Nusselt num- 
ber variation in parallel -plate channel 
with uniform wall heat flux. Specular 
reflection coefficient, 1; ratio of specific 
heats, 1.4; Prandtl number, 0.73. 



Ratio of slip to average velocity, u $ /u 


Figure 4. - Value of definite integral for parallel-plate 
channel for any value of slip to average velocity ratio. 
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Figure 5. - Variation of dimensionless wall-to- 
bulk temperature difference at heated section 
entrance for slip flow in parallel -plate channel 
with uniform wall heat flux. 
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(c) Temperature-jump coefficient, 0.4. 


Figure 6. - Wall temperature ratio in thermal entrance 
region for flow in para I lei -plate channel with uniform 
wail heat flux and different values of temperature-jump 
coefficient. Symmetry number, 2. 
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(c) Temperature-jump coefficient, 0.4. 


Figure 7. - Wall temperature ratio in thermal entrance 
region for flow in parallel-plate channel with uniform 
wall heat flux and different values of temperature-jump 
coefficient. Symmetry number, 1. 
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Rarefaction Rarefaction 

parameter, parameter, 



Dimensionless axial distance, fx/2L)/RePr 

(a) Symmetry number, 2. <b) Symmetry number, 1. 

Figure 8. - Wall temperature ratio in thermal entrance region for flow in a parallel-plate channel with uniform wall heat 
flux. Specular reflection coefficient, 1; ratio of specific heats, 1.4; Prandtl number, 0.73. 



Dimensionless axial distance, <x/2L)/RePr 

(a) Symmetry number, 2. (b) Symmetry number, 1. 

Figure 9. - Variation of Nusselt number along parallel -plate channel for uniform wall heat flux and different values of 
temperature-jump coefficient. 
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Dimensionless axial distance, tx/2L)/RePr 


(a) Symmetry number, 2. 


(b) Symmetry number, 1. 


Figure 10. - Variation of Nusselt number along parallel-plate channel for uniform wall heat flux. Specular reflection 
coefficient, 1; ratio of specific heats, 1.4; Prandtl number, 0.73. 



Figure 11. - Physical model and coordinate system for round tube. 
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Ratio of slip to average velocity, u s /u 

Figure 12 - Values of definite integral for any value of 
slip to average velocity ratio in round tube. 
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Dimensionless axial distance, (x/r )/RePr 
o 

Figure 13. - Variation of Nusselt number along round 
tube for uniform wall heat flux and different values 
of temperature-jump coefficient. 
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Figure 14. - Variation of Nusselt number along round 
tube for uniform wall heat flux. Specular reflection 
coefficient, 1 : ratio of specific heats, 1.4; Prandtl 
number, 0.73. 



Figure 15. - Wall temperature ratio in thermal entrance 
region of round tube for uniform wall heat flux and 
temperature-jump coefficient of zero. 



